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nanospheres, within a framework that incorporates the spatial variation of the optical fields within 
the sphere. Our recent theory of exchange dipole spin wave modes of the sphere provides us with 
eigenvectors. When properly normalized, these eigenvectors allow calculation of the absolute cross 
section of various modes which contribute to BLS spectrum. We then present explicit calculation 
of the BLS spectrum associated with the first few dipole/exchange spin wave modes with emphasis 
on their relative intensity. 
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I. INTRODUCTION 

As lithographic techniques have rapidly developed in 
the last decade, it has become possible to assemble or- 
dered magnetic elements on the nanometer scale. These 
have attracted interest as potential candidates in mag- 
netic data-storage devices. Spin wave excitations control 
the dynamic response of the magnetization in the linear 
response regime, and the speed of real devices at least for 
small amplitude motions of the magnetization. This pro- 
vides motivation for theoretical and experimental studies 
of spin waves in magnetic wires, dots and other nanoscale 
structurea^i^i^. 

In our previous paper—, we provided an analyt- 
ical method to calculate eigenvalues and eigenvec- 
tors of dipole-exchange spin waves in ferromagnetic 
nanospheres, and we explored the properties of these 
modes. A question is how one may access these modes 
experimentally. We argue here that Brillouin light scat- 
tering spectroscopy (BLS) offers the means of probing the 
dipole-exchange spin waves of very small spheres. In this 
paper, we develop the theory of BLS from a single metal- 
lic nanosphere, and illustrate with explicit calculations 
that one may access the spectrum of dipole/exchange 
spin waves modes by this means. 

As we know, there are two complementary experimen- 
tal techniques to study spin waves in small structures 
such as nanoscale particles or spheres: ferromagnetic res- 
onance (FMR) spectroscopy and Brillouin light scatter- 
ing (BLS). FMR uses microwaves to excite spin waves, 
but the microwave wavelength is much larger than the 
radii of nanospheres. It is the case as well that their 
radius is far smaller than the skin depth of microwaves, 
which is typically a micron for materials of current in- 
terest. Thus, FMR produces a uniform exciting field 
inside nanospheres. This makes it impossible to excite 
the full spectrum of dipole-exchange modes of such an 
object. Only the uniform mode is excited in this circum- 
stance. However, the laser photon used in BLS has a 
small skin depth at optical frequencies. This is typically 
10-20 nm in the ferromagnetic metals of interest. This 
is comparable with the radius of nanoscale spheres, and 
the resulting non-uniform optical field can thus, in prin- 



ciple, excite a full spectrum of dipole-exchange modes. 
It is then of interest to develop an explicit description 
of the scattering cross sections associated with the vari- 
ous modes, and to obtain insight into the factors which 
control these cross sections. This is the purpose of the 
present paper. We shall see below from the calculations 
presented in the present paper that can indeed excite and 
study higher order dipole exchange modes through the 
BLS method. We shall conclude that BLS is a powerful 
tool which one can provide information about the nature 
of the spin wave spectrum of nanoscale spheres. This will 
be true of samples of less simple shape as well, of course. 
A quantitative theory of the BLS spectrum of more com- 
plex objects is quite non trivial, unfortunately. While 
our calculations apply only to small spheres, our prin- 
cipal conclusion that the method can excite the dipole- 
exchange spin wave spectrum of nanoscale magnetic en- 
tities is quite general, in our view. 

From the quantum point of view, the process of BLS 
in ferromagnets is viewed as an interaction between 
light photons and magnons. From the classical point of 
view, the thermal fluctuations of magnetization inside 
the sphere will change the dielectric tensor of the fer- 
romagnetic material in a time dependent manner. The 
inelastic scattering event has its origin in the scattering 
of the laser light from these dynamic fluctuations in the 
dielectric constant. We use a semi-classic method here 
to address this issue. 

In this paper, we present the theory of BLS from spin 
waves in ferromagnetic nanospheres, and we calculate the 
cross section of BLS from several low lying dipole ex- 
change modes. We shall see that for small spheres, the 
cross section for exciting particular higher order dipole- 
exchange modes is appreciable, indeed larger than that 
of the uniform mode for sufficiently small spheres. 



II. THE STRUCTURE OF THE THEORY 

We consider a metallic ferromagnetic nanosphere with 
radius in the range of lOnm to 60nm, with vacuum out- 
side. Since the linear dimensions of the sphere are large 
compared to the lattice constant, this allows us to use 
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continuum theory to describe the spin motions in the 
sphere. The ferromagnetic nanosphere is magnetized 
uniformly along z-direction with saturation magnetiza- 
tion Ms- At small but finite temperatures, the ther- 
mally excited magnetization can be written by M(r, t) = 
MsZ -\- rn{r, f), where m(r, t) is a small amplitude distur- 
bance associated with thermal spin waves. As noted ear- 
lier, the eigenvectors which describe the dipole-exchange 
spin waves of the sphere can be obtained from our pre- 
vious paper—. These allow us to describe the spin wave 
contribution to rn{f, t) provided that the eigenvectors are 
properly normalized. We remark that the prescription 
for normalizing such eigenvectors has been discussed only 
quite recentlj^.In the regime where dipolar contributions 
to the excitation energy are appreciable, the prescription 
for normalizing the eigenvectors is not so obvious. 

The fluctuation of magnetization causes a time- 
dependent change of the dielectric tensor of the mate- 
rial, which we write as 5e^^{f, t) =^ K^v\m\{r, t). We 

A 

retain terms only first order in the transverse magneti- 
zation, to generate a description of BLS by processes in 
which a single spin wave quantum is created or destroyed. 
When applied to static magnetization arrangements, the 
term just described controls Faraday rotation. An early 
theoretical analysis of the spin wave BLS spectra of 3d 
ferromagnetic films shows that a fully quantitative ac- 
count of surface and bulk spin wave contributions follow 
by choosing X^i^a = Ke^^x, where e^j^a is the Levi Civita 
tensor—, appropriate to bulk materials of cubic crystal 
structure. We adopt this picture here, so that we have 



5e^y{r,t) = K 
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xmx{r,t). 



(1) 



As just discussed, thermal fluctuations produce time 
and space varying components to the magnetization, 
which in turn lead to fluctuations of the dielectric ten- 
sor. Now we turn to the description of the light scattering 
of electromagnetic radiation by dielectric tensor fluctu- 
ations in the sphere. We describe the electromagnetic 
fields in the vicinity of the sphere by Maxwell's equa- 
tions: 



VxE = — — , 
c ot 



(2) 



where uj is the frequency of the incident laser. Upon 
inserting this into Maxwell's equations, we have 



V X E{r,t) = ie-B{r,t), 
c 



(6) 



V X B{f,t) ^ ^ie-E{r,t) + —,pff{?^t). (7) 
c c 

We have defined J'Jf{r,t) = ^ ^ 5e^,(f, i). 

Upon using the relation between the dielectric tensor and 
the magnetization as in Eq. ([T]) , we obtain the form of the 
effective current density, given by 

J'J^ (^^' = ii E -^M->^"^A (r, t)E, (f, t) . (8) 

Since the thermal fluctuations are very small in ampli- 
tude, we will treat the effect of the term involving J^^^ in 
Eq.(l7|) in the Born approximation. Thus, we may replace 



E^,{r,t) by Elr{T',t) in Eq.®, with E), ' the amplitude 
of the incident field. Thus, we have 



AmA(f,t)i?(o)(f,t). (9) 



The problem is now that of computing the radiation 
fields generated by the effective current density within 
the sphere. Of interest to us are radiation fields in the 
far zone, as r ^ oo. We may proceed by following the 
approach outlined by Jackson^. Far from the sphere, the 
radiation fields may be described as a superposition of 
fields of TM and TE character. In the far zone, the TM 
component of the radiation field has the magnetic field 
perpendicular to f, the unit vector in the radial direction, 
while the electric field has both a radial and transverse 
component. The TE mode has the character of E and B 
interchanged. Given the radial component of E in the far 
zone, one can derive expressions for the transverse com- 
ponents of both E and B for the TM mode. Furthermore, 
one can derive equations satisfied by r- E everywhere: 



Attuj 



r-.r^f,{r < R), (10a) 



V X B 



I dD 
c~dt' 



(3) 



where within the sphere, 

with e is the complex dielectric constant of the mate- 
rial from which the sphere is fabricated. Since (5e^^(r, t) 
varies slowly in time with respect to the laser field, from 
Eq.(|3]) we have, 

^^^^^ = -ieujE^(T, t) - tuJ2Se,ir, t)E,{r, t), (5) 



[V^ + ^]ir-E)^0,{r> R). 



(10b) 



Similar statements can be made regarding the TE fields. 
One can derive equations satisfied by f- B, and given B 
in the far zone, one can derive the remaining components 
of B. The combination r • B satisfies: 



[V' + ^e{L.W-B) 



L- ,{r < R), (11a) 



](r.i?) = 0,(r>E), 



(lib) 
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where L = ir x V is the orbital angular-momentum op- 
erator of quantum mechanics. 

We can calculate the radial components of the electro- 
magnetic fields by using Green's functions. By having 
solutions for the radial components of B and we can 
generate a complete description of the TM and TE mul- 
tipole mode fields in the far zone, where the fields have 



radiative character. The analysis proceeds very much 
along the lines given in Jackson's text^, so we only pro- 
vide the final results. The transverse field components in 
the far zone are of interest, since these control the radial 
flow of outgoing scattered energy. The expressions for 
these field components in the far zone may be written 



= E (^)i7(h:T) [^r^^rf -)l,ii,„.(^, ^) - Ar)rf„-)x,y,™(^, ^)], (i2a) 
E^' = E (^)II(}tT) [^'""'^'-^^^ ^) + Ar)r[--)L,y,™(^, ^)], (i2b) 

- E MWTT) [^'""'^'™'^^^'™^'' ^) + ^r"^r[--)L,y,.(^, ^)], (i2c) 



where — lo/c^ and we note the angular momentum 
operators — and Le = -^^-§^- The quantities 

'^Un^'^ ^'^'^ '^fm^ '"^'^ viewed as the (Z, m) components 
of the effective current densities, with origin in the ther- 
mal fluctuations of the magnetization, as we have seen. 
The source rl"^^-* generates scattered radiation of TM 

character, while r[^^ generates scattered radiation of 
TE character. We nave the forms 

rS'^ = / '^Ji{kr)YCmiO, V)r- ^"{r), (13a) 

and 

rl^"" = / d'rMkr)L-.rff{f)Y*^{e,^)^ (13b) 
where k = fco\/e is the (complex) wave vector of elec- 



tromagnetic radiation inside the sphere. The imaginary 
part of k controls the penetration depth of the incident 
laser radiation. 

In Eq.(12), the factors A^^*"^' and Af^^ have the 
following physical interpretation. The effective cur- 
rents just described create scattered radiation inside the 
nanosphere; this radiation must be transmitted through 
the surface of the sphere to the outside. These quantities 
are the transmission coefficients for TM and TE waves of 
partial wave character I, respectively. We need to discuss 
properties of these transmission coefficients to appreciate 
the reason why the numerical results in section III dis- 
play the systematics we shall see. We begin by quoting 
the general expressions for a|"^*'''* and A^^'^'', and then 
we pause to explore their behavior in the limit of interest 
to us. 

We express the factors as 



,(TM) 4nkok\R^ ^ ji{kR)hf\{kR) - ji.,{kR)h'i'\kR) 

A; = ( ), (14a) 

c €HiikoR)jiikR) - hY\koR)Gi{kR) 

^ .(TE) Ank^R.ji{kR)h['_\{kR)-ji^^{kR)h['\kR) 

and A = ( — ). (14b) 

c Hi {koR)ji (kR) - h\^> {koR)Gi (kR) 



We have introduced Gi{x) — xj/_i(x) — lji{x) and Hi{x) = xh\]}j^{x) — lh'i^\x), where ji{x) is the spher 



ical Bessel function of the first kind, and /ij (x) is the 
spherical Hankel function of the first kind. The wave 
vector k in Eq.(14) is again the (complex) wave vector of 
electromagnetic radiation within the material from which 
the nanosphere is fabricated. Thus, the magnitude of the 
product kR is roughly the radius of the nanosphere di- 
vided by the optical skin depth S. As noted in section 
I, this ratio will be in the range of unity, for samples of 



interest. The quantity fco is the free space wave vector of 
electromagnetic radiation of frequency co and for spheres 
whose diameter is near the smaller end of the the range of 
10-60 nm we shall have koR ^ 1. Thus, of interest is the 
relative magnitudes of a['^*^'' and a['^'^'' when koR ^ 1. 
Upon using small argument forms for the relevant Hankel 
functions, one finds in this limit 



.(TM) ^ {koR}'+^ ^ ink^R ^ MkR)h\'],ikR)-ji^,ikR)h['\kR) ^^ 
' "(2/-l)!!^ c ^ ljiikR) + GiikR)/e ^ ' 

, ,iTE) ^ (fcofi)'+^ ^ 'i7^k^RJlikR)h'^\{kR)-jl^^{kR)h['\kR) ^^ 
=(27^^^( lMkR)+G,ikR) 



In Eq.(15a) and Eq.(15b), the factors in curly brackets 
are close in magnitude to each other. We thus see that 
in the limit that the wavelength of the laser radiation 
is large compared to the radius of the sphere, the in- 
tensity of TE radiation with angular momentum I ( its 

strength is controlled byA^"'"'^-' ) will be stronger than the 
TM radiation with angular momentum I ( its strength is 
controlled byA|"^*^'') by the large factor of {ki)R)~^ . This 

follows since r|"^^^ and r["^^ are roughly the same order 
of magnitude. This observation will be very important 
when we discuss the relative scattering intensities of the 
various dipole-exchange spin wave modes below, in the 
limit the radius of the sphere is small. 

We have one final step. In the definition of the ef- 
fective current density J'^-f^{r,t) given in Eq.([9]), we re- 
quire an expression for the incident laser field E^^^^ (r) 
inside the nanosphere. If we consider a plane wave in- 
cident on a conducting sphere of arbitrary radius, the 
expression for this field is cumbersome. Here we ex- 
ploit the fact that our sphere has radius R small com- 
pared to the wavelength of the laser light in vacuum. In 
this limit, to excellent approximation, we may assume 
that the sphere is placed in the spatially uniform field 
i^C) = zE'^^^exp{—iujt). In the region just outside the 
sphere, we are far into the near zone, to use language ap- 
propriate to radiation theory, so the spatial form in the 
field outside the sphere may be taken to be that given by 
the equations of electrostatics. However, the radius of the 
sphere is not small compared to the optical skin depth, so 
we must solve for the electric field inside through use of 
the full electromagnetic theory. We surely have V • i? = 



inside the sphere, so the electric field obeys the Helmholtz 
equation 



V2^(f) + k^E{r) = 



(16) 



everywhere inside the sphere, where as above k^ = 
(uj/c)'^e with e the optical frequency dielectric constant. 
Outside the sphere, in response to the incident field we 
have the classical electrostatic field of a polarized sphere. 
In the near zone and this has components 



E: 



2a 



El 



■ sin 0, 



whereas inside the sphere the fields have the form 



(17a) 



b d 
r dr 



(rji(fcr)) sin( 



(17b) 



Application of the boundary conditions is straightfor- 
ward, and yields 



6=?- 



2en{kR) + f^{rn{kR))\i 



(18) 



The expressions for optical fields inside a conducting 
nanosphere under conditions similar to those considered 
here will prove useful for a variety of applications. Thus 
we pause by arranging the results into a more useful for- 
mat. We define 



R 



e + 2 



r {e-l)ji{kR) + sin{kR) 



(ji (kr) + sin(fcr)) 



(19a) 
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and A/ = 



R 



r {e- l)ji{kR) + sm{kR) 
I 



(3ji(fcr) — sin(fcr')). 



(19b) 



Of course, inside the sphere, the field is invariant under 
rotation around the z axis. The two non zero components 
of the field in cylindrical coordinates are given by 



(/ + A/(r) cos(20)) 



(20a) 



and Ep{r) 



A/(r)sin(26i). (20b) 



In the limit that the skin depth is large compared to the 
radius of the sphere, the function f(r) approaches unity 
everywhere, and A/(r) approaches zero. We are then 
left with the elementary expression for the fields inside a 
dielectric sphere exposed to a spatially uniform electric 
field. When the skin depth is comparable to or smaller 
than the radius of the sphere, the use of Eqs.(20) allows 
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us to account for the spatially non uniform character of 
the laser field inside the sphere. It should be remarked 
that in the numerical calculations we report below, the 
magnetization of the sphere is always taken parallel to 
the z axis, and we will explore the Brillouin cross section 
when the laser field is not necessarily parallel to the z 
axis. Of course, in this case one may apply a simple 
coordinate rotation to the field components in Eqs.(20) 
to generate expressions for the exciting field inside the 
nanosphere. 

We now have all the ingredients in place to obtain an 
expression for the Brillouin scattering cross section. One 
evaluates the Poynting vector in the far zone, integrates 
it over solid angle, and divides the result by the energy 
per unit time which illuminates the sphere. One finds the 
following for the total Brillouin scattering cross section: 



ctbls 



V ^ (\\iTM)-^(TM)f2 



(21) 



The expression in Eq.(21), written out in full, is a 
quadratic form in the transverse magnetization com- 
ponents, with terms proportional to the combination 
m^(f,t)mv{r',t) integrated over the volume of the 
sphere; one has an integral over both r and f*. One 
takes a thermal average over this quadratic form. The re- 
sulting expression provides one with the Brillouin cross 
section, integrated over all solid angle, and integrated 
also over all spin wave modes which contribute to the 
Brillouin signal. One can break the correlation functions 
< mf^{f,t)m,y{r',t) > down into sums over each of the 
spin wave eigenmodcs of the sphere, and thus extract 
from Eq.(21) the Brillouin cross section integrated over 
solid angle for each individual mode. When we refer to 
the Brillouin cross section of an individual mode, we in- 
clude both the Stokes and the anti Stokes component. 
The algebra associated with the decomposition just de- 
scribed is complex, and since the resulting forms are of 
little general interest, we omit the details here. We refer 
the reader to the discussion given in ref. [6] for a descrip- 
tion of how this is done. It is the case, as noted in this 
paper, that in the regime where dipolar contributions to 
the spin wave excitation energy are of comparable mag- 
nitude, the algorithm for normalizing the eigenvectors is 
not so obvious, as noted above. The means of doing so 
are derived in ref. [6] . 

We turn next to our numerical studies of the Brillouin 



scattering intensities of the various normal modes. 



III. RESULTS AND DISCUSSION 

We first begin with some general comments. It is the 
case that the nature of the spin wave modes in ferro- 
magnetic spheres and related sample shapes is a classic 
topic in magnetism. Motivation for the early theoreti- 
cal studies was provided by the first generation of FMR 
experiments, which were largely carried out on spheri- 
cal samples of yttrium iron garnet (YIG). These samples 
had linear dimensions that were macroscopic. As a con- 
sequence, the long wavelength spin wave modes excited 
by the microwave fields used in FMR could be described 
by magnetostatic theory with exchange ignored. In a 
remarkable paper that is now classic. Walker presented 
an analytic theory of the magnetostatic spin wave modes 
of uniformly magnetized ferromagnetic ellipsoids^. The 
special case of the sphere was discussed by Fletcher and 
Bell, who provide a very large and interesting array of 
useful formulae for descrbing the eigenfrequencies and 
eigenvectors as welli^. As the radius of the sphere is 
made small, to the point where the radius lies in the few 
nanometer range, exchange enters importantly in the de- 
scription of the spin wave eigenmodes. This is particu- 
larly the case for the materials of current interest, which 
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are the 3d ferromagnets or alloys fabricated from these 
elements. By virtue of the itinerant character of their 
ferromagnetism, along with their high Curie tempera- 
tures, the exchange stiffness D is very large. The first 
description of the exchange/dipole spin wave modes of 
a ferromagnetic sphere has been developed only recently 
by Arias and the present authors^; this paper extends the 
analysis of ref.[9] and ref.[10] to include exchange as well 
as dipolar interactions between the spins. A complete 
discussion of the response of such a sphere to spatially 
inhomogeneous microwave frequency magnetic fields is 
given as well. 

As discussed in ref. [1], a measure of the importance 
of exchange in the description of the spin wave modes 
of a ferromagnetic sphere is the dimensionless param- 
eter Tex = (^)^' where R is the radius of the 
sphere. Consultation of Fig. 4 of ref.[l] shows that the 
contribution of exchange to the excitation energy of spin 
waves in spheres is substantial even for values of r^x as 
small as 0.2. The calculations presented below are for 
Fe nanosphcres, for which AirMg = 2 x W^gauss, and 
D ~ 2.5 X 10^^ gauss ■ cm?. For these parameters, and 
for an Fe sphere with radius of 10 nrn, we have r^x = 1- 
We call the reader's attention to the Appendix of ref.[l] 
where simple analytic expressions are provided for the 
excitation energies of spin waves in the limit of strong 
exchange, with the first correction from dipolar interac- 
tions included. As noted in ref. [1], this expression ac- 
counts nicely for the excitation energies calculated from 
the full theory for a rather wide range of r^x ■ 

In this paper, we refer to the dipole exchange spin wave 
modes by the quantum number scheme used in ref.[l], 
which is adapted from ref. [9] and ref. [10]. First, suppose 
we consider the niagnetostatic limit. Then for each eigen- 
mode of the ferromagnetic sphere, the magnetic dipole 
field generated outside the sphere by the spin motion has 
an angular dependence controlled by a single spherical 
harmonic 11, m(^, v)- Inside the sphere, the magnetic po- 
tential which describes these fields is not simple, if wishes 
to express it in terms of spherical harmonics. It is the 
case that each mode can be labeled by the two quan- 
tum numbers Z, m associated with the spherical harmonic 
just mentioned, and a radial quantum number n which 
may be viewed as a radial quantum number, to make an 
analogy with quantum mechanics. Hence, in the mag- 
netostatic limit, each mode is labeled by the three num- 
bers {l,m,n). The uniform mode, which is the mode 
seen in FMR when the sphere is exposed to a spatially 
uniform microwave field, is the (1,1,0) mode in this no- 
tation. When exchange is added, it is no longer true 
that the magnetic potential outside the sphere is pro- 
portional to a single spherical harmonic. However, as 
the limit that the exchange stiffness D is allowed to ap- 
proach zero, each mode of the exchange/dipole spectrum 
smoothly and continuously reduces to a particular mag- 
netostatic mode which can be labeled by the scheme just 
described. Thus, we label each dipole/exchange mode 
of the sphere by the same three quantum numbers em- 
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FIG. 1: For the case where the electric field of the laser is 
parallel to the magnetization, we plot the cross section ratios 
o"2io/o"iio and craio/fno as a function of the radius of an Fe 
nanosphere. The influence of exchange, as a function of ra- 
dius, for these and other dipole/exchange modes is given in 
Fig.4 of ref.[l]. 




FIG. 2: We show the dependence of the cross section ratios 
o"2io/o"iio and cr20o/o"iio on the angle 9i between the magneti- 
zation of the sphere and the electric fleld of the incident laser 
beam. The calculations are for an Fe sphere whose radius is 
10 nm. 



ployed for the mode in the limit D ^ Q. We recall ref.[l] 
that with no spin pinning at the surface of the sphere 
(the case considered here), the FMR mode is unaffected 
by exchange, and remains a uniform mode even with ex- 
change present. We now turn to our results. 

We have carried out calculations of the BLS cross sec- 
tions for the (110) mode, the (2m0) modes and the (3m0) 
modes, and also the (410) mode. Since there is little in- 
terest in the absolute BLS cross section for a single, iso- 
lated sphere, we shall present the result in the form of 
cross section ratios of the form cr;mo/fiio- Our principle 
conclusions are summarized in Fig.l and Fig. 2, for the 
case where the incident photons have the energy of 2.5 
eV. In Fig.l, for the case where the electric field vector 
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of the incident radiation is aligned with the magnetiza- 
tion of the sphere, we show two of the cross section ratios 
as a function of the radius of the nanosphere. We have 
taken material parameters appropriate to Fe, but clearly 
the trends are robust, and not affected by details of the 
material parameters. We see that as the radius of the 
sphere is reduced, the scattering power of the (210) mode 
increases dramatically relative to the FMR mode, while 
that of the (310) mode remains roughly similar to that 
of the FMR mode. 

What is responsible for the behavior of the cross sec- 
tions in Fig.l is the operation of what we may call a 
quasi selection rule. First, for the FMR mode, one may 
show that all of the coefficients r[^^ vanish identically, 
and the only non zero contribution to the cross section 
comes from ■ Thus, we can say that in the limit 

when fco-R <C 1, the dominant contribution to the BLS 
cross section for the FMR mode (that from F^^^') is 
"silent", and the FMR mode scatters only by virtue of 
the contribution r^^*^\ which as we see from Eqs.(15) 
leads to a contribution to the cross section smaller than 
that from the leading term by the factor of (fcoi?)^. For a 
nanosphere with a radius of 10 nm, we have koR ~ 0.13, 
whereas when we reach 50 nm. radii, kgR is sufficiently 
large that the quasi selection rule is not operative. For 
the (210) mode, T^^^ is substantial for the smaller 
spheres. The quasi selection again suppresses the cross 
section of the (310) mode in the limit of small radii. 

In Table I, we tabulate the cross section ratio for sev- 



eral modes, and for sphere radii of 10 nm, and also 50 
nm. Wo sec that for the larger sphere, the (110) has 
the largest cross section, but the (310) mode is nearly as 
strong. As the sphere is reduced in size the quasi selec- 
tion rule dominates, and the I = 2 modes dominate the 
spectrum, to have very large cross sections save for the 
m = 2 case, where the cross section is very small; we 
insert a zero into the Table for this reason. 

The angular dependence of the cross section ratios is 
also of interest. In Fig. 2, we show the cross section ratios 
as a function of the angle between the incident laser field 
and the magnetization for the two modes which domi- 
nate the spectrum for small spheres. We see that in the 
limit koR <C 1 the large cross section for the (210) mode 
remains large for all angles, whereas the scattering from 
the (200) mode is suppressed as electric field of the inci- 
dent laser is rotated from parallel to the magnetization, 
to the equatorial plane of the sphere. Experimental ob- 
servation of this angular dependence would be of great 
interest, as a test of the picture put forward in this paper. 

Insight into the origin of the quasi selection rule which 
controls the results above in the limit of small radii can 
be obtained from the limiting forms displayed in Eqs.(15) 

and in the structure of the quantities r|^^ and r["^^ 
defined in Eqs.(13). For the form of the coupling be- 
tween the light and the spin system described earlier, 

and for the case where the incident electric field is aligned 
with the magnetization of the sphere, one has the explicit 
forms: 



(22a) 



and t\^J^^ =iAj^ ^rji{kr)f{r)Yi^{e, sm{e)m^{f). 



(22b) 



In these expressions, the integration is over the volume of 
the nanosphere, we have defined A = \^ujKE'^^'' /A-K{e + 
2)] , and the components of the transverse magnetization 
associated with a given spin wave mode are expressed in 
cylindrical coordinates. 

From the structure of the eigenvectors given in ref. [1] , 
one may see that each of the dipole/exchange spin wave 
modes of the sphere has well defined parity, even or 
odd, under reflection through the equatorial plane of 
the sphere. The parity of the mode may be written as 
(— 1)'+™. Now the spherical harmonic which appears in 
Eq.(22b) along with the factor of sin(0) also have parity 
(— 1)'+™ under this reflection, whereas the parity of the 
integrand in Eq.(22a) is (—1)'+™+^. Then we also note 
that for a given mode with label {lo,mo,n) the angular 



integration over ip renders all the quantities in Eqs.(22) 
zero except for the contributions for which m = toq. 

Thus, we can see from the previous paragraph that 
the term which provides the leading contribution to the 

cross section, F^"^^-* vanishes for the (110) mode, which 
is the mode excited in FMR. It is the case that for this 
mode, the only non zero contribution to these coefficients 
is Til ■ The FMR mode is thus forbidden to scatter 
in lowest order, in the limit fco-R <C 1. When we turn 
to the I = 2 modes, for the (200) mode F^q^^ is non 
zero and substantial, and for the (210) mode the same is 
true of r^^'^''. For the (220) mode, the azimuthal inte- 
grations render all the leading terms, {T^i".^^^}, zero and 
we obtain our flrst substantial contribution from fI^ . 
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Cross Section Ratio 


R — lOnm 


R — 50nm 


Cross Section Ratio 


R — lOnm 


R — 50nm 


at 6/ = 






at 9i = n/2 






C"2(Jo/ <^ 111) 


29.20 


0.09 


O"20o/o"ll0 


0.09 


0.001 


C21o/f'"lll) 


32.43 


0.03 


O"21o/o"ll0 


32.43 


0.02 


O'22o/f"110 





0.01 


O"22o/o"110 





0.005 


csoo/fiio 


1.13 


0.49 


O"30o/o"110 








O"31o/o"llo 


0.37 


0.07 


O"310/o"ll0 


1.51 


1.70 


O'32o/f"lll) 


0.01 


0.001 


O"32o/o"110 


0.007 


0.001 


f 330 /fill) 








O"33o/o"110 








O'41o/o" 110 





0.02 


O"41o/o"110 





0.01 



TABLE I: For spheres with radii of 10 nm and 50 nm, respectively, we tabulate the cross section ratios aimo/'^iio for the 
various dipole/exchange modes considered. We provide results for the case where the incident laser field is parallel to the 
magnetization {6i — 0) and the case where it is perpendicular {9i = 7r/2). Zeros are inserted into the table for modes with 
very small excitation cross sections. 



The cross section for this mode is thus very small, in 
the Umit kgR < 1. For the (310) mode, r^'['^^ is forbid- 
den by reflection symmetry through the equatorial plane, 
and the first non zero contribution is from F^"^*^-*. The 
cross section for this mode is thus comparable to that for 
the (110) mode, the FMR mode. We can continue on 
to understand the small cross sections for the (320) and 
the (330) mode. By the time we reach the (410) mode, 
the angular momentum content of the spin wave is suffi- 
ciently high that the scattering cross section is small. 

From the calculations presented in this section, we can 
see that BLS spectroscopy can allow one to access se- 
lected higher order dipole exchange modes, whereas fer- 
romagnetic resonance spectroscopy probes only the single 
FMR mode in the limit of small sphere radii. 

IV. CONCLUDING REMARKS 

From the calculations presented here, we can see 
that BLS spectroscopy of the spin wave modes of small 
nanosphcres can provide one with access to certain higher 
order dipole/exchange modes. In the limit of small sphere 
radii, in the sense described above, we see that the I = 2 
modes are predicted to dominate the spectrum. 

The analysis presented here explores the BLS cross sec- 
tion for the various spin wave modes, for a single, isolated 
nanosphere. Of course, any experimental study will nec- 
essarily explore an array of such objects. Indeed, Stein- 
muUer et al. have reported BLS studies of Fe nanosphcres 
deposited on a GaAs surfaceii. It is difficult to comment 
in detail on their results, since the film consisted of clus- 
ters of Fe spheres arranged in a disordered fashion on 
the surface. We can comment on one issue, however. 
These authors observe Stokes/anti Stokes asymmetries 
in their spectra, and argue this provides evidence that 
the spin excitations they study should be regarded as 
collective excitations of the array of spheres, which inter- 
act through dipolar coupling. We recall that Stokes/anti 
Stokes asymmetries in BLS spectra taken from the sur- 



faces of ferromagnetic crystals and in films were reported 
many years ago, and theory based on the physical pic- 
ture used here provides a fully quantitative account of the 
data^^'^'^'-^"'. A question raised with these experiments is 
whether the Stokes/anti Stokes asymmetry is truly ev- 
idence for collective mode behavior, or whether in fact 
such asymmetries are found for the scattering of light 
from spin waves in a single, isolated sphere such as that 
studied here. In section III, we presented only total cross 
sections, integrated over both the Stokes and anti Stokes 
features in the spectrum. We remark that in the course 
of this investigation, we have inquired into the possible 
occurrence of Stokes/anti Stokes asymmetries for the sin- 
gle sphere, to find the answer in the negative. Thus, we 
reinforce the interpretation offered in ref.[ll]. 

It would be highly desirable to see BLS studies from 
ordered lattices of ferromagnetic nanosphcres, as opposed 
to the disordered films studies in ref.[ll]. Then, of course, 
one would observe collective spin wave modes of the lat- 
tice of spheres where the intersphere interaction is dipolar 
in nature, if one assumes that neighboring spheres do not 
touch. We remark that we have developed a theoretical 
description of such collective modes, for the case where 
both exchange and intrasphere dipolar interactions are 
comparable in strengtbi^. In our paper, one finds explicit 
calculations which describe the nature of such modes. It 
will be a straightforward matter to adapt the description 
of the BLS spectrum of an isolated sphere set forth in 
the present paper to the description of scattering from 
the collective modes. Experimental data on ordered ar- 
rays of spheres will provide motivation for this extension. 
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